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Wc present a magneto-optical study of the three-dimensional topological insulator, strained HgTe 
using a technique which capitalizes on advantages of time-domain spectroscopy to amplify the signal 
from the surface states. This measurement delivers valuable and precise information regarding the 
surface state dispersion within <1 meV of the Fermi level. The technique is highly suitable for the 
pursuit of the topological magnetoelectric eflfect and axion electrodynamics. 



Spin 1/2 particles exhibit a counterintuitive property 
that their wavefunction acquires a tt phase upon 360° ro- 
tation. If spin and orbital degrees of freedom are mixed 
in a particular way, the momenta of electrons in a crys- 
talline lattice feel important effects of this tt 'Berry's 
phase', which can lead to a new phase of matter whose 
description requires a fundamental redress of the theory 
of semiconductors [1-3], and probably many other mate- 
rials classes [4-6] . These topological insulators exhibit an 
odd number of metallic surfac;c bands with helical spin 
texture [7-9] surrounding the nominally insulating bulk 
and display characteristic suppression of backscattering 
from step edges and nonmagnetic impurities [10, 11]. We 
present magneto-optical measurements deep in the ter- 
ahertz frequency regime exploring the charge dynamics 
of surface states in high-mobility strained films of HgTe 
[12-15]. Using a time-domain technique, we detect a 
strong magneto-optical signal which is dominated by sur- 
face bands. This information reveals precise details of the 
low-energy excitations and momentum-energy dispersion 
of the helical metallic surface state. 

The observation of the quantum spin Hall effect 
(QSHE) in CdTe-HgTe-CdTe quantum wells [16, 17] rep- 
resents a significant advance in the ability to robustly 
segregate electronic currents of opposite spin, an effect 
which paves the way to new applications for spintronics 
and fault-tolorant quantum computation. In bulk, HgTe 
possesses the band inversion property due to spin-orbit 
interaction, a prerequisite condition for the QSHE, but 
also needed to realize a throe dimensional topological in- 
sulator phase. Unfortunately, the Fermi level appears 
directly at the intersection of two bands, rendering the 
bulk scmimctallic and leaving the surface states and the 
expected topological aspects of this material obscured by 
the low-energy bulk excitations. However, when HgTe 
is compressivcly strained against a CdTc substrate, the 
bulk band intersection becomes fully gapped in response 
to the lowered symmetry, permitting isolated access to 
helical surface bands [12 1-5]. Recently, strong evidence 
for the existence of the n Berry's phase, surface bands. 



and associated zero Landau level were observed as a 
quantum Hall effect of surface states, whose existence 
were subsequently verified using angle-resolved photoe- 
mission [15]. 

In order to probe directly the surface states of strained 
HgTe, we have developed a method of acquiring time- 
domain terahertz magneto-optical data using a home- 
built superconducting magnet in a flow cryostat, il- 
lustrated in Figure la. Complementary to frequency- 
domain measurements performed on the same sample at 
higher temperature (>50 K) [18], which were interpreted 
in terms of a thermally-activated carriers in the bulk, 
our measurements were performed at low temperature, 
where the terahertz response is due to a combination of 
intrinsically doped bulk carriers and topological surface 
bands [15]. At 4.35 K, the magnetic field and incident 
polarizer angle 6p were set, and one direct pulse plus 
one echo, due to internal reflection inside the substrate, 
wore collected for analyzer angles 6a subtending 400° in 
10° steps. The pulse and echo were then partitioned into 
Ei{t, 9a, B) and E2{t, 9a, B) as shown in Figure lb, and 
separately Fourier transformed to obtain the electric field 
amplitudes Ei{uj, 9a, B) and E2{lo, 0a, B). 

When the pulse passes through or reflects from the 
film, the electric field direction can rotate due to off- 
diagonal elements of the dielectric tensor or possibly 
the predicted topological magneto-electric effect [19]. 
Accounting for the horizontally-polarizing detector and 
emitter antenna of the spectrometer, the electric field 
amplitude in pulse i at the detector is 

Ei = Eg cos 9p cos 9a cos{9a — 9p + 9F,i) (1) 

where Eg is the electric field generated by the source 
and 9F,i is the complex, sample- induced Faraday rotation 
of the ith pulse. Figure ld,e shows an example fit of 
equation (1) to the measured intensity \E2(i^,9a,B)\'^ of 
the second pulse for fields of ±1.4 T. The change in shape 
of the intensity distribution around 35 cm~^ is caused 
by a change in the complex Faraday angle Op^2 at this 
frequency, which we attribute to an absorptive transition 
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FIG. 1: (Color online) Time-domain measurement of the 
normal-incidence Kerr angle, (a) Illustration of the magneto- 
optical THz apparatus. 6p was set to 51.2° throughout, (b) 
Subset of the angle-dependent measurements of the 1st and 
2nd pulses. Vertical lines show the cuts used to separate the 
pulses into Ei and E2. (c) Measured (red dashes) and 

the fit from eqn (1) (blue lines) used to determine the Fara- 
day angle ^^,2(1^) for +1.4T and -1.4T. (d) Pictorial rationale 
for the expression 0f,2 — 0f,i = Ok- (e) The measured 6k for 
several magnetic field values. 



among surface state Landau levels. The fitting procedure 
was applied to both pulses to uniquely determine 9p^i and 
9f,2 at each frequency and field value measured. 

The simultaneous acquisition of multiple pulses with 
separate histories of contact with the film provides an 
advantage in accurately determining the polarization ro- 
tation due to the film. Any polarizing effects of the 
polyethylene cryostat windows, as well as inaccuracy in 
repositioning the polarizer following an angle sweep is 
exactly the same for the two pulses, so taking the differ- 
ence in the two Faraday angles 0k=9f,2 — Op,! cancels 
these extraneous effects [23]. A separate measurement of 
a reference CdTe substrate showed negligible rotation in 
magnetic field, and this normal-incidence Kerr angle 6^, 
shown in Figure le and 2a-b is precisely the rotation of 
polarization induced when the pulse reflects at normal 
incidence from the substrate side of the film (Fig. Id), 
and the accuracy to which it can be determined is greater 
than for the Faraday angles separately. 

The field-induced changes in Ok are typical of cy- 
clotron resonance (CR) behavior, for which the dynami- 
cal conductivity appropriate for photon angular momen- 
tum either parallel (+) or antiparallel (-) the momentum 



of the incident photon is 



(2) 



where o;^ is the cyclotron frequency, the Drude weight, 
and 7^ is the inverse lifetime of the jth band. The el- 
ementary response functions in Equation (2) are then 
plugged into the expression 



tan 6ii 



iZons{(T+ - (T-) 



(l + Zoa+)(l + Zoa_) 



to fit the measured Kerr angle. Here, is the measured 
substrate index of refraction and Zq is the impedance 
of free space (see Supplemental Materials [24]). While 
we explore the possibility of a multi-component CR be- 
low, the resonance at each field is fit and satisfactorily 
described by a single effective uJc, 7, and ujf. This pro- 
cedure results in the parameters summarized in Figures 
2c-e. The sign of the rotation indicates that the carri- 
ers involved in the CR are electron-like, and fits to the 
field dependence give fi'^ ~ 0.9 meV, a field-linear cy- 
clotron frequency with huc/B = 2.92 ±0.02 meV/T, and 
total Drude weight Tiojf — 79.0 ± 1.7 meV. Using the 
relation fiB = uJc/li and the nearly field-independent 7 
(Figure 2d), these fits imply a very high carrier mobility 
/i — 34, 220 cm^/V- s, consistent with previous transport 
[15] and optical [18] measurements. 

In the limit of small carrier concentration, which as we 
will see below is the relevant limit in the context of this 
discussion, the 2D Fermi-surfaces are isotropic, and wj- 
and at low fields are uniquely determined by the Fermi 
momentum, kp, and the group velocity at the Fermi- 
energy, Vp (see Supplemental Materials) 

■ eB v'p 
uji — — ^ . 

Knowledge of wj- and uj^ therefore permits determination 

of the parameters Vp and kp, which in turn provides 
valuable and precise information on the band structure 
at and around the Fermi energy. 

In our 70 nm films, the conduction band is quantized 
and well-separated in energy due to confinement in the z 
direction, and at low temperature any possible bulk con- 
tributions come from a small number of two-dimensional 
electron pockets. Unlike bulk, the surface contributions 
to the Drude weight are present for all values of the chem- 
ical potential, due to their gapless nature. Because the 
surface bands are at higher filling {kp > kp) and are 
more steeply dispersing {vp > Vp) than the conduction 
bands, these states not only always contribute to the CR, 
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FIG. 2: (Color online) Properties of the surface-state cy- 
clotron resonance. Real (a) and imaginary (b) Kerr angle 
6k{ijJ,B) determined from terahertz ellipsometry. (c) The 
experimentally determined ujc (red squares) and a linear fit 
(purple lines) showing the B-linear CR behavior, ujc cor- 
responds to the inflections in Re&x(i^), as demonstrated by 
the experimental frequency derivative KedOx /du, shown as a 
false color plot in the background. (d,e) Fit parameters j and 
ijf. The error bars arise from scatter in the determined CR 
parameters. 



but also always contribute more strongly to the magneto- 
optical signal than a set of bulk states with the same 
Fermi surface area. 

Before we address the implications of possible contri- 
butions from the bulk, we start with the simplest as- 
sumption, namely that the only contributions to the 
observed Kerr rotation originate from the two surfaces 
of the film with approximately the same charge car- 
rier concentration. This is motivated by the magneto- 
transport showing a quantized Hall effect that results 
from the 2D Dirac-like topological surface states with 
densities 4.8xl0"cm-2 and 3.7xl0"cm-2 for the CdTe 
and vacuum interfaces (respectively) and negligible con- 
tribution from the bulk [15]. Since with this assumption 
hujj = 39.5 ± 0.9 meV for the two surfaces, combina- 

5x10^ m/s and fcf, = 0.201 



tion with uic/ B gives Vp —b.'i 
nm~^ (see Fig 3a). The corresponding carrier concentra- 
tion per surface is n2D = kp/i^T^) = 3.2 x lO^^cm"^ in 
good agreement with the high field transport data [15] 
and strongly evidencing a surface-dominated origin of 
the CR signal. The remaining difference is within the 
observed variations from one cool-down to another of 
the same sample, which we attribute to molecular ad- 
sorption at the surface of the film at cryogenic tempera- 



tures. Comparing these kp values with the quantum well 
bandstructure in [15], the level of the chemical potential 
should be positioned above the conduction band bottom 
with 25% of the Drude weight arising from the bulk states 
of the film, with the remainder due to the surface bands. 
Attributing only 75% of the observed Drude weight to the 
surface states would result in: Vp =5.09x10^ m/s and 
kp = 0.174 nm~^, as indicated by the green dotted hues 
in Fig. 3a. Angle-resolved photoelectron spectroscopy 
shows a linear dispersion down to 1 eV below the Fermi 
energy with a velocity 4.3x10^ m/s [15]. The somewhat 
higher value obtained from the optical data is expected, 
since this probes the velocity at the Fermi energy which 
is closer to the light conduction bands hybridizing with 
the 2D surface states. 

A more detailed perspective of the surface-dominated 
cyclotron resonance can be attained through comparison 
of these results to a surface state model dispersion rela- 
tion = Ep)p+Ak+C2k^ [20]. This form accommodates 
the topologically protected Dirac point (DP) through the 
A term, but allows for significant deviations from the 
ideal conical dispersion through the C2 term. By vary- 
ing these parameters, one can interpolate continuously 
between a pure Dirac cone (C2 — ^ 0) and a parabola 
{A — >■ 0) emanating from the Dirac point at zone cen- 
ter and energy Ep)p lying inside the bulk gap. All of 
the curves of this type with the same values oi kp and 
Wi? fit equally well to our data, as shown in Figure 3b. 
However, in order to place the DP inside the bulk gap, 
estimated to be ~ 26 meV at zone center [15], the A 
parameter must be very small, A < 50 meV-nm, to be 
consistent with our data. The resultant parameter C2 is 
therefore quite large in comparison to other topological 
materials [7-9, 20], indicating a relatively rapid departure 
from the ideal conical dispersion near the DP in strained 
HgTe. The large second-order term leads one to suspect 
that higher order isotropic terms {k^, fc^, ...) may be 
necessary for detailed analysis of certain experiments. 

While we do not measure the Fermi energy directly, 
the simplified dispersion analysis above puts the chemical 
potential within 40 meV of the Dirac point. This differ- 
ence can be overcome through gating in an appropriate 
experimental design, and the chemical potential can be 
tuned near the Zeeman split zero Landau level. This 
would fulfill the conditions required to observe the pre- 
dicted topological magnetoelectric coupling effect [19, 21] 
as a quantized Kerr rotation. This long-sought effect 
bears close mathematical analogy to high-energy parti- 
cle theory, permitting one to use terahertz spectroscopy 
to study the properties of an 'axion' domain wall [22], 
and the methods developed here are highly suited for 
this pursuit. 

We have studied the low energy electrodynamics of 
topological surface states of strained HgTe using a 
novel time-domain magneto-optical spectroscopic tech- 
nique providing Kerr angle spectra to very low energies 
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FIG. 3: (Color online) Dispersion of the surface state band, 
(a) Constraints on Fermi parameters kp and vp consistent 
with the magneto-optical data. Black line shows the con- 
straint from the observed Uc/B, and purple dashes indicate 
the error bars on this quantity. Creen solid curve is a con- 
straint from the observed cj/ assuming the CR is due solely to 
two identical surface states at the surface and interface, with 
solid orange curves indicating the error bars. Dashed orange 
and green curves show the effect of breaking this assumption 
and allowing 25% of the Drude weight to be attributed to the 
confined conduction band states, (b) Possible realizations of a 
minimal surface state model in [20] which are consistent with 
our observations for different choices of A and C2. The Fermi 
level fiF, shown in black, intersects the surface and confined 
bulk conduction subbands. (c) Dispersion of the surface state 
quasiparticles determined in this work with 75% surface state 
signal, A — 25 meV-nm, and C2 = 890 meV-nm^. 



(<1 meV). The method allowed us to obtain the param- 
eters describing the topological surface states near the 
Fermi energy, until now not resolved by other experi- 
mental techniques, namely free carrier spectral weight, 
quasi-particle scattering rate, cyclotron frequency, Fermi- 
velocity and Fermi-momentum. Taken together with the 
requirement that the Dirac point must lie inside the bulk 
gap, our results imply that the surface bands of strained 
HgTe are markedly nonconical. 
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the National Center of Competence in Research (NCCR) 
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Materials and methods 



High mobility, 70 nm HgTe film sample was grown on a thick CdTe substrate using molec- 
ular beam epitaxy, and characterized by X-ray, ARPES, and d. c. transport measurements. 
The finite area of the film, 4x5 mm, limits the lowest measurable frequency in our measure- 
ments. Time-domain terahertz traces were collected with a commercial system employing a 
mode-locked laser and photoconducting antennae. 



Magneto-optical THz measurements 

Equation (1) of the main text is an example of Malus' law, and is straightforward to 
derive. 

For a given angular momentum state of the photon (-|- or -), the transmission coefficient 

is 



for the first pulse and 



for the second pulse, where 



2 

= (3) 

are the usual Presnel coefficients for transmission and refiection at the substrate/ vacuum 
interface, 0^ is the complex phase shift acquired through traversal across the substrate and 
s, /, and V stand for substrate, film, and vacuum, respectively. 
The coefficients which include the thin film are 

^'^""'^ ^ n, + l + Zoa^ 

where Zq = ^/Jio/eo ~ 376.63... is the impedance of free space. The Faraday angle of the 
ith pulse 9F,i is 

tan%^-^ ^^'+~^^'- (7) 

ti,+ + ^j,- 



3.30 




FIG. 1: Index of refraction of a CdTe substrate measured using terahertz spectroscopy and used 
to fit the normal incidence Kerr angle 9k to the cyclotron resonance formulae (Eqn 2, main text). 
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FIG. 2: Fit quality of the CR formula to the measured Re^i^ for positive fields. 

Using a trigonometric identity, 

tan Ok = tan {eF,2 - Of,i) = ^ ^^or^s(cT+ a ) ^ 

This relation connects the model conductivity to the measured normal-incidence Kerr angle. 
The frequency-dependent index of refraction of the substrate was determined from separate 
measurements and is shown in Figure 1. The quality of fit used to deduce Uc, ujf, and 7 is 
shown for positive fields in Figure 2. 
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Drude weight and cyclotron frequency, general considerations 



For isotropic planar dispersion, the Drude weight of a band i is 

The frequency u of electrons executing closed orbits of /c-space area A satisfies [1] 

1 27re de , 
B h^ojdA ^ ' 

where the cyclotron frequency uic corresponds to the first harmonic {n — 1). 

In our case of an isotropic 2D Fermi surface, 

A = T:kl (11) 



We use the fact that de/dk = Tiv is the group velocity, so that 

— = 27r— 
de hvp 



= 2.-^ (12) 



from which 



hojc = eB^. (13) 



Surface state model 

The surface state model is given in reference [2] as 

es{k) = Edp + Ak + C2k\ (14) 

The Fermi wavevector is 

A;^ = -^+VA' + ^C,{^ip-Ej,p) ^ ^^^^ 
2C2 

and the Fermi velocity is 

nv'p = A + 2C2k'p. (16) 
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